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Abstract: In this paper, the prediction and investigation on drag coefficient of non-spherical particles is presented 


applying artificial neural network. The performance between BP (Backpropagation) model and RBF (Radical Basis 


Function) model is compared. The RBF model is employed to predict the drag coefficient of non-spherical particles 


with higher efficiency and less error. The simulation results are compared with the experimental results in the 


literature by using RBF model. It reveals that artificial neural network can be applied to the prediction on the drag 


coefficient of non-spherical particles. 
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0 引言 


固体 颗粒 在 连续 性 流体 中 的 输 运 是 


= NA 


普遍 现 


象 。 针对 固体 颗粒 的 输 运 过 程 , 研究 者 们 从 实验 和 数 


值 模拟 两 方面 i 


非 球形 颗粒 的 存在 。 
在 ， 其 运动 特性 与 球 


行 了 大 量 的 研究 03。 但 是 ， 研 究 多 
集中 于 球形 颗粒 , 而 实际 工业 过 程 中 不 可 避免 的 会 有 
由 于 非 球形 颗粒 形状 因 
颗粒 有 着 较 大 的 差异 。 
研究 非 球形 颗粒 在 连续 怕 


子 的 存 
因此 ， 


流体 中 的 运动 状态 十 分 必 


要 。 上 颗粒 在 流体 的 输 运 过 程 中 , 气 固 之 间 动 量 交换 对 
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颗粒 运动 行为 的 控制 起 着 至 关 重 要 的 作用 。 研 究 者 们 
通过 构建 蝶 力 模型 的 方法 来 描述 气 固 之 间 的 动量 交 
换 , 其 中 电力 系数 的 选择 和 计算 对 电力 计算 的 准确 与 
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否 有 着 不 可 忽视 的 作用 。 

根据 目前 的 研究 结果 可 以 发 现 , 针对 颗粒 在 流体 
自由 运动 时 上 忠 力 系数 的 研究 , 实验 方法 是 一 种 比较 
容易 的 方法 ,研究 者 们 主要 采用 沉降 管 和 风 洞 实验 等 
方法 对 非 球 形 颗粒 的 点 力 系 数 进行 预测 3。 但 是 ， 
采用 实验 测量 的 方法 也 存在 不 足 , 例如 采用 风 洞 方法 
时 存在 实验 空间 较 大 , 调试 安装 工作 较 复杂 ， 人 力 物 
力 耗 费 巨 大 等 。 

近 些 年 来 ， 人 工 神经 网 络 模型 取得 了 长 足 的 发 
展 , 由 于 其 具有 能 够 精确 估算 复杂 问题 ， 比 唯 算 模型 
效率 高 , 对 于 不 完整 的 有 品 输 入 数据 具有 更 高 的 效率 
等 优点 ， 因 此， 人 工 神 经 网 络 模型 正 广泛 应 用 于 不 同 
学 科 的 研究 之 中 ， 可 以 有 效 避 免 弥补 实验 方法 的 不 
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究 方向 为 多 相 流 流动 。 


足 。 目 前 ， 在 流 态 化 领域 


相关 研究 已 有 相关 报道 。 研 究 者 们 采用 
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经 网 络 模型 进行 


预测 和 分 析 流 化 床 内 氮 氧 


神经 网 络 模型 


攻 物 的 生成 过 程 , 泄露 现象 
的 检查 和 预测 以 及 固体 垃圾 气 化 等 过 程 , 并 


取得 了 民 


好 的 效果 [9。 但 采用 人 工 神 经 网 络 模型 对 气 固 暇 力 


系数 进行 研究 的 工作 ， 还 鲜 有 报道 。 
综 上 所 述 , 采用 神经 网 络 模型 预测 非 球形 颗粒 的 
气 固 忠 力 系 数 的 研究 还 很 少见 。 因 此， 本 文采 用 了 神 


经 网 络 模 型 预测 非 球 


颗粒 的 气 固 忠 力 系数 。 首 先 ， 
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接着 ， 采 用 RBF 模型 对 非 球形 颗粒 


型 的 预测 效果 ， 


1 BP 人 工 神经 网 络 模型 示意 图 
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Fig. 1 Schematic diagram of BP neutral network model 


BP 算法 的 提出 在 一 定 程度 上 解决 了 多 


数 训练 难 的 问题 ， 但 是 其 自身 也 有 局 


BP 算法 需要 较 多 参数 ， 如 网 络 层 数 、 


经 网 络 模 型 的 预测 误差 ; 


气 固 忠 力 系 数 进 


行 预测 和 对 比分 析 , 给 出 了 不 同 


力 系数 随 雷诺 数 的 变化 情况 。 
1 人 工 神 经 网 络 模型 


人 工 


脑 结构 及 其 功能 的 信 
种 计算 工具 ， 由 若干 简单 、 
研究 
输入 信号 的 动力 系统 响应 ， 


与 过 程 信息 。no 


1.1 BP 人 工 神经 网 络 方法 


息 ,处 理 系统 9 神 


定义 的 神经 元 组 成 。 


即 ANN), 通常 简称 为 神经 网 络 模型 ， 
神经 网 络 的 启示 下 建立 的 数据 处 理 模 型 ， 
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子 情况 下 ,由 


经 网 络 模型 (Artificial Neural Network， 


是 一 种 在 生物 
则 在 模仿 人 
经 网 络 模型 是 一 


高 度 内 在 连接 的 、 由 许多 


神经 网 络 
能 够 学 习 


方法 通过 对 外 在 
高 度 非 线性 关系 


1974 年 ，Werbo 等 人 (9 首次 提出 了 用 于 神经 网 


络 学 习 的 BP 算法 ， 为 多 层 神经 网 络 的 学 习 训练 与 实 
现 提 供 了 一 种 切实 可 行 的 解决 途径 。1986 年 ， 


Rumelhart 和 McCelland 等 
差 反 向 传播 算法 进行 了 详尽 的 分 析 , 解决 了 多 属 


科学 家 


网 络 的 学 习 问题 ,进一步 反 
网 络 的 拓扑 结构 包括 输入 层 、 隐 含 层 和 输出 层 ， 能 够 


Ez 了 BP 


对 多 层 网 络 的 误 
经 


BP 


算法 的 发 展 。 


在 不 知道 输入 输出 其 体 数 学 表达 式 的 情况 下 , 通过 学 


习 来 存储 复杂 的 映射 关系 。 在 BP 网 络 中 ， 


参数 的 学 


习 通 常 采用 误差 反 向 传播 算法 , 数据 自 输 入 层 经 隐 含 
层 逐 层 向 后 传播 , 网 络 的 连接 权 值 通过 输出 层 经 过 中 
间 各 层 逐 层 向 前 , 沿 着 误差 性 能 函数 梯度 的 反方 向 修 


正 , 原理 采用 最 速 下 降 法 。BP 网 络 模 型 如 图 


1 所 示 。 


元 个 数 及 权 值 ， 
BP 算法 采 | 


限 性 。 


展 网 络 参 
首先 ， 


层 的 神经 


每 一 


且 缺 乏 有 效 的 参数 选择 方法 ， 其 次 ， 
j 最 速 梯度 下 降 的 优化 思想 ， 在 解决 实际 


问题 时 , 会 经 常 陷入 到 局 部 极 小 值 中 ， 导致 算法 陷入 
局 部 最 优 ， 再 次 ，BP 算法 对 于 样本 的 依赖 性 较 高 ， 


如 果 样 本 集合 代表 性 差 、 存 在 大 量 矛 盾 样 本 和 元 余 样 


本 ， 网 络 的 性 能 就 会 较 差 。 最 后 ， 对 于 一 些 复杂 网 络 


的 优化 问题 ， 
几 个 小 时 。 


BP 算法 


1. 2 RBF 基 神 经 网 络 方法 


因此 ， 随 着 神 将 网 络 模型 的 发 
们 逐渐 提出 了 新 的 神经 网 络 模型 以 弥补 


1988 年 ， 基 于 生物 神 


经 元 的 局 部 响应 特 


受到 学 习 速 率 的 限制 需要 花费 


展 ， 研 究 者 
其 不 足 。 
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生 ， 


Broomhead 和 Lowe 将 径 向 基 函 数 引 入 到 了 神经 网 络 


的 设计 中 中 


， 提 出 了 径 向 基 神 


经 网 络 RBF 算法 。 后 


来 ，Jackson 和 Park0415 分 别 于 1989 年 和 1991 年 对 


RBF 算法 在 非 线性 连续 函数 


了 论证 。RBF 神经 网 络 是 一 种 3 
基本 工作 原理 是 : 利用 RBF 构成 的 非 线 性 隐 合 
间 对 低 维 的 输入 矢量 进行 投影 , 将 数据 变换 到 高 维 
间 中 去 , 以 使 原来 线性 不 可 分 的 问题 能 够 变 得 线性 可 
分 , 最 终 得 到 线性 的 输出 层 。 


的 基本 结构 示意 图 。 
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2 RBF 人 工 神经 网 络 模型 示意 


上 的 一 致 逼近 性 能 进行 
层 的 前 向 网 络 ， 
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图 2 为 径 向 基 神 经 网 络 
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Fig. 2 Schematic diagram of RBF neutral network model 


RBF 网 络 有 很 快 的 学 习 收 敛 速度 ， 一 个 很 重要 
的 原因 在 于 其 结构 简单 ， 有 较 快 的 收敛 速度 , 没有 学 
习 隐 含 层 权 值 的 过 程 , 省 去 了 误差 在 网 络 中 逐 层 传递 
的 耗 时 过 程 。RBF 网 络 的 研究 与 应 用 ， 标 志 着 神经 
网 络 真正 走向 实用 化 过 程 , 目前 已 被 广泛 应 用 于 非 线 
性 函数 逼近 、 模 式 分 类 、 控 制 系统 建 模 、 时 变数 据 分 
析 和 故障 分 析 诊 断 等 工程 领域 。 

本 文 的 输入 参数 为 雷诺 数 Re 和 形状 因子 P， 预 
期 输出 值 为 气 固 电力 系数 。 采 用 已 有 文献 中 的 部 分 实 
验 值 作 为 学 习 和 训练 数据 , 预测 不 同形 状 因子 下 气 固 
忠 力 系数 随 雷诺 数 的 变化 情况 。 

2 结果 与 讨论 

2.1 BP 与 RBF 算法 计算 结果 对 比 

本 节 首 先 应 用 BP 和 RBF 两 种 方法 对 文献 中 不 
同形 状 因子 的 抽 力 系数 的 预测 结果 进行 了 对 比 。 
Pettyjohn 和 Christiansen 等 人 09 的 实验 分 别 给 出 了 形 
状 因子 为 0.670、0.806、0.846、0.906 和 1.000 时 上 归 
力 系数 分 布 的 实验 结果 。 本 文 将 形状 因子 为 0.670、 
入 0.846 和 1.000 时 的 实验 结果 用 于 学 习 和 训练 ， 通 过 
一 、 ”训练 , 预测 以 上 5 种 形状 因子 下 遇 力 系数 的 分 布 ， 
RJ] 与 实验 结果 进行 对 比 。 对 比 结果 如 图 3 所 示 。 当 形状 
因子 为 0.906 时 ,无 论 是 采用 BP 方法 还 是 RBF 方法， 
预测 得 到 的 忠 力 系数 与 实验 结果 均 较 为 接近 , 但 采用 
BP 方法 要 比 采 用 RBF 方法 消耗 的 时 间 较 长 ， 如 表 1 
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3 BP 与 RBF 方法 预测 结果 对 比 
Fig. 3 Comparison of simulation results by BP and RBF models 


图 4 给 出 了 采用 不 同形 状 因子 颗粒 时 , 预测 结果 


.三 所 示 。 当 形状 因子 为 0.806 时 ， 如 图 3 (b) 所 示 , 采 
-一 用 RBF 方法 的 预测 结果 与 实验 结果 符合 更 好 。 
采用 BP 方法 和 RBF 方法 对 不 同形 状 因子 时 电 
力 系数 预测 的 计算 时 间 对 比如 表 1 所 示 。 从 表 中 可 以 
明显 看 出 ， 在 相同 运行 环境 下 ， 采 用 BP 方法 的 计算 
时 间 平 均 在 80s 左右 ， 采 用 RBF 方法 的 计算 时 间 在 
8s 左右 ， 是 BP 方法 的 10%。 由 此 可 见 ， 采 用 相同 数 
据 量 作为 学 习 和 训练 数据 预测 相同 形状 因子 情况 下 
的 虑 力 系 数 时 ，RBF 方法 的 计算 效率 更 高 。 
表 1BP 和 RBF 计算 时 间 对 比 


Table 1 Comparison of computing time 


esl 


形状 因子 1.000 0.906 0.846 0.806 0.670 


与 实验 结果 的 平均 误差 分 布 情况 。 从 图 中 可 以 明显 看 
出 ， 采 用 RBF 方法 的 预测 误差 始终 小 于 10%， 预 测 
误差 分 布 较为 稳定 。 而 采用 BP 方法 预测 结果 误差 分 
布 波动 较 大 ， 特 别 是 当 形 状 因子 为 0.670 时 ，BP 方 
法 的 预测 误差 接近 21%， 预 测 误差 分 布 介 于 10% 和 
20% 之 间 ， 误 差 较 大 ， 预 测 结果 的 可 信 性 较 低 。 


综 上 所 述 ， 采 用 


BREF 方法 预测 结果 更 接近 于 实 


验 值 ， 预 测 误 差 小 ， 


效率 高 ， 对 于 不 同形 状 因子 颗粒 


BP 83s 87s 85s 84s 84s 
RBF 8.41s 8.93s 9.01s 8.90s 8.14s 


气 固 忠 力 系数 的 预测 是 一 种 较为 合理 且 有 效 的 方法 ， 
J 以 为 不 同形 状 因子 情况 下 的 点 力 系数 的 选取 提供 
更 合理 的 理论 依据 。 因 此 ， 本 文 将 采用 RBF 的 方法 
对 不 同形 状 因子 颗粒 气 固 上 忠 力 系数 进行 预测 和 分 析 。 
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图 4 BP 与 RBF 预测 误差 分 析 
Fig. 4 Error analysis of BP and RBF models 

2. 2 不 同形 状 因子 电力 系数 预测 

由 上 文 对 比 可 知 ，RBE 方法 可 以 很 好 的 预测 不 
同形 状 因子 情况 下 的 气 固 点 力 系 数 ， 并 且 耗 时 较 短 。 
Pettyjohn 和 Christiansen 等 人 的 实验 结果 共 提 供 了 5 
组 实验 数据 , 本 文选 取 其 中 3 组 用 于 训练 和 学 习 ,， 3 
预测 用 于 训练 和 学 习 的 形状 因子 情况 下 的 电力 系数 ， 
与 实验 结果 进行 了 对 比 , 对 比 结果 如 图 5 所 示 。 从 图 
中 可 以 发 现 , 当 形 状 因子 为 1 时 预测 结果 与 实验 结果 
对 比较 好 , 说 明 对 用 于 训练 和 学 习 的 数据 进行 预测 可 
以 得 到 理想 的 预测 结 

由 图 可 见 , 电力 系数 随 着 雷诺 数 的 增加 而 逐渐 降 
低 ， 当 雷诺 数 增 加 到 100 左右 时 ， 暇 力 系数 值 基本 保 
持 不 变 ， 但 是 当 雷 诺 数 达 1000 时 ， 暇 力 系数 略 有 增 
加 。 预 测 结果 与 实验 结果 平均 误差 为 8%， 对 比 结 
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图 5 RBF 网 络 预测 忠 力 系数 与 文献 实验 数据 比较 (学 习 数 据 ) 


Fig. 3 Comparisons between RBF simulation results and 


experimental data (test data) 
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为 了 进一步 说 明 RBF 预测 方法 的 合理 性 ， 本 文 
将 训练 学 习 到 的 模型 对 研究 工 况 中 的 另外 2 组 形状 
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因子 下 忠 力 系数 进行 了 预测 , 并 与 实验 结果 进行 了 对 
比 ， 对 比 结果 如 图 6 所 示 。 从 图 中 可 以 明显 发 现 , 预 
测 结 果 与 实验 结果 吻合 较 好 。 在 雷诺 数 转折 位 置 
CRe=100) 处 ， 预 测 结果 与 实验 结果 略 有 偏差 ， 但 


整体 预测 效果 与 实验 结果 较为 接近 , 平均 预测 误差 约 
为 7% 和 89%。 
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图 6 RBF 网 络 预测 忠 力 系数 与 文献 实验 数据 比较 (预测 数据 ) 


Fig. 6 Comparisons between RBF simulation results and 


experimental data (predicted data) 


综 上 所 述 ，RBF 预测 结果 较为 准确 ， 可 以 为 不 

同形 状 因 子 情况 下 的 忠 力 系数 进行 预测 , 并 为 非 球形 
颗粒 忠 力 系数 的 计算 提供 了 合理 的 理论 依据 。 
图 7 给 出 了 采用 RBF 模型 对 不 同形 状 因子 颗粒 
气 固 电 力 系数 进行 预测 的 结果 。 文 中 分 别 预 测 了 形状 
因子 为 0.5、0.6、0.7、0.8、0.9 及 1.0 时 电力 系数 的 
分 布 情况 。 从 图 中 可 以 看 出 ， 在 不 同形 状 因子 下 ， 熏 
力 系 数 随 雷 诺 数 的 变化 情况 相似 : Re<10 时 ， 上 忠 力 系 
数 与 雷诺 数 基本 呈 线 性 变化 , 随 着 雷诺 数 的 增加 蝶 力 
系数 逐渐 减 小 ， 当 10<Re<100 时 ， 电 力 系数 变化 趋 
势 逐 渐 趋 于 平缓 ; 当 Re>100 时 ， 电 力 系 数 基本 保持 
不 变 。 对 比 不 同形 状 因子 下 电力 系数 的 分 布 可 以 发 
现 ， 当 Re<1 时 ， 电 力 系数 彼此 较为 接近 ， 当 Re>1 
时 , 不 同形 状 因子 下 , 遇 力 系数 的 差距 逐渐 变 得 明显 ， 
并 且 随 着 雷诺 数 的 增加 形状 因子 对 昌 力 系数 的 影响 
逐渐 增强 。 当 Re>100 时 ， 不 同形 状 因子 下 电力 系数 
随 雷诺 数 的 变化 逐渐 减 小 ， 数值 上 略 有 增加 , 但 基本 
保持 不 变 。 
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7 RBF 神经 网 络 方法 预测 点 力 系数 随 形状 因子 的 变化 情况 
Fig. 7 Influence of shape factor on drag coefficient by RBF 


model 
3 结论 


本 文采 用 人 工 神经 网 络 模型 , 对 非 球形 颗粒 气 固 
电力 系数 进行 了 预测 及 分 析 , 并 将 模拟 结果 同文 献 中 
实验 数据 进行 对 比分 析 。 结果 表明 ， 人 工 神 经 网 络 可 
以 用 于 非 球 形 颗粒 气 固 中 力 系 数 的 预测 研究 ,主要 结 
论 如 下 : 

(1) BP 模型 和 RBF 模型 预测 结果 与 实验 结果 
对 比 发 现 ， 采 用 RBF 模型 的 预测 结果 与 实验 结果 更 
加 接近 ， 误 差 较 小 ， 计 算 效 率 较 高 ; 

(2) 采 用 RBF 方法 预测 不 同形 状 因子 下 电力 系 
数 的 分 布 , 可 以 发 现 , 不 同形 状 因子 下 电力 系数 随 雷 
诺 数 的 变化 趋势 基本 相同 ; 

(3) 在 相同 雷诺 数 下 ， 随 着 形状 因子 的 减 小 外 
力 系数 逐渐 增 大 。 

综 上 所 述 , 采用 神经 网 络 的 方法 可 以 较 好 的 预测 
不 同形 状 因 子 下 的 气 固 上 忠 力 系数 , 同时 预测 结果 表明 
形状 因子 对 忠 力 系数 的 影响 不 可 被 忽略 , 本 文 的 工作 
对 非 球形 颗粒 气 固 上 忠 力 模型 的 预测 和 构建 具有 一 定 
的 参考 价值 。 
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